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Abstract: The feasibility of producing superheavy nuclei in proton evaporation channels was systematically studied 


within the dinuclear system (DNS) model. Due to the Z=114 proton-shell, one can synthesize Fl isotopes in proton 


evaporation channels. We only considered the case of evaporating one proton first and then n neutrons in this work, 


other cases were ignored due to the small cross-section. The production cross sections of unknown isotopes 


290,291 Fy 


in 38S+255Es reaction are the highest compared with ??Ti42^Np and ?^Cr-2??Pa reactions, and the maximum cross 


sections are 1.1 and 15.1 pb, respectively. ?S-2* Es is a promising candidate to approach the island of stability as the 


radioactive beam facilities are upgraded in the future, and the production cross sections of ”!~?**F] in that reaction are 


estimated to be 3.2, 6.0, 4.0, and 0.1 pb, respectively. 
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Introduction 


The synthesis of superheavy nuclei has always been an 
important task and challenge in nuclear physics, which is sig- 
nificant to expand the nuclear map, investigate the origin of 
heavy elements, and test the shell model. At present, the re- 
search on the synthesis of superheavy nuclei can be divided 
into two directions. One is to synthesize new elements to ex- 
plore the charge limit of superheavy nuclei [13]. The other 
is to move towards to the double magic nuclei FI, which 
is the center of the island of stability ^". Because the £ sta- 
bility line bends towards the neutron axis, the currently pro- 
duced superheavy nuclei based on fusion-evaporation reac- 
tions are neutron-deficient, and very far away from the island 
of stability. That requires us to strive towards the develop- 
ment of radioactive beam facilities and the search for new 
production mechanisms. 


Several models have been developed to study the fu- 
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sion mechanism and predict the production cross-sections 
of superheavy nuclei. The semi-classical models, such as 
the multidimensional Langevin-type dynamical equations !?l, 
the nuclear collectivization model [9] the fusion-by-diffusion 
modell!’ and the dinuclear system modell-?!, have been 
successfully applied to calculate the evaporation-residue 
(ER) cross sections of superheavy nuclei. The microscopic 
models like TDHF approach!!*5! and ImQMD model !!6-171 
can effectively describe the dynamical dissipation during the 


fusion process. 


In this work, we attempted to produce superheavy nuclei 
in proton evaporation channels. To this end, we have devel- 
oped the de-excitation part in the DNS model and found the 
production of Fl isotopes in proton evaporation channels is 
considerable. 


1 The DNS model 


Within the framework of DNS, the ER cross sections of 


superheavy nuclei at the incident energy E, m. are calculated 


c.m. 
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QJ + LT(E. mJ) 
A =È 


x Pon(Eom.> J) Wr(Ecnm,, J). (1) 


Here, T, Pox, and Wr are the transmission, fusion, and sur- 


Ogg CE, ss ) = 


vival probability, respectively. 

It is difficult to obtain an analytical expression for the 
transmission probability due to the complexity of the ac- 
tual Coulomb barrier, and it can be obtained by using the 
WKB method after making a parabolic approximation of the 
Coulomb barrier”; 


T(Eems J) = jio 


1 


cm B= 2 zo + DII 


d B. 


1 十 exp{ 一 LE 


a 
(2) 
Here, fic is the width of the parabolic barrier and Rp de- 
fines the position of the barrier. Considering the multi- 
dimensional character of the realistic barrier, the barrier dis- 
tribution function f(B) should be introduced, which is taken 
as the asymmetric Gaussian form!”!. 
The interaction potential after considering deformation 


can be written as 


V(R, fi. f». 0;, 0) = Vo(R, fi. f». 1, 02) 
1 
+V (R, fj, By, 91,92) + 5€ 


where Vo and Vy are Coulomb potential and nuclear poten- 
tial, given by Wong's formula and double folding method, re- 
spectively. Bis is the static deformation (for nuclei | and 2), 
usually taken as the quadrupole deformation. fi; is the dy- 
namical deformation. Assuming that the deformation energy 
is proportional to the mass number, i.e. C, f?/C, f? = A,/A;, 
only one deformation parameter f = fj, + f; can be used 
to represent the dynamical deformation. C,» are the coeffi- 
cients that characterize nuclear hardness and given by 


3 Ze 

=(A-1) (A+DRvo- oa QA) : (4) 

Here, Ry is the nuclear radius, 4 is the deformation number, 

taking 2 for quadrupole deformation. o is the surface tension 

coefficient and given by 4z R?c = a, 423, where a, is taken 
as 18.32 MeV. 

The distribution probability P(Z,, N,.1) for the frag- 

ment (Z,,N,) at time ¢ is determined by solving the master 


- BY + ZCP- SY. 
3) 
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Here, dz x, discribes the microscopic dimension corre- 
sponding to the macroscopic state (Z1, Nj). Aq; and Ag, 
are the quasifission and fission rates, calculated by the one- 
l Ot) = Ve*lais the local 


temperature, obtained from the Fermi gas model, where e* 


dimensional Kramers equation”! 


is the local excitation energy of the dinuclear system, and a 
is the level density parameter. Wz Ni ZN (Wz, NZ Nn’) is 
the transition probability from state (Z,, Ni) to (Z,, Nj) (or 
from (Z4, N,) to (Z,, N1)), which can be written as: 


Ten Z1 NE, Z, , Ni, Ei;D) 
hr? dz n dz! ,Ni 


x È |< Zp Ni, Epi VOIZ, Ny, Ei > [. (6) 


sah 
H 


W; NZ, Ny (D = 


Here, i represents the remaining quantum numbers. E, de- 


notes the local excitation energy. T 


‘mem 1S the memory time 
[22] 


and given by 


(Aj, Ej, Aj, Ej, th= (7) 
inay : MeeVee < a RK”? 


The compound nucleus formed by fusion reaction has 
high excitation energy and deexcites by emission of y-rays, 
light particles and fission. According to Weisskopf’s evapo- 
ration theory), the evaporation width of particle v can be 


written as: 
D(E',J) = Qsj41)————— 
ESI) = Qs Do EU 
E*-B,-6-6,-+ 
x ent" — B, — ó, — €) (Ede, 
0 
(8) 
where B,, m,, and s, are the binding energy, mass, and spin 


ofthe particle, respectively. The pairing correction 6 is set to 
be 12/ VA, 0, and 一 12/ VA for even-even, odd-A, and odd- 
odd nuclei, respectively. 6, is the neutron correction energy. 
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Fig. 1 (color online) The ratios B/B, of nuclei 2867297 NI, 
286-297 g| 285-26Mc, and 78°-?°7Ly are denoted by the cir- 


cles, squares, triangles, and pentagrams, respectively. 


If the neutron number of the compound nucleus is odd, 9, = 
12/1 A, and 6, —0 in other cases. p is the energy level density 
and calculated by the Fermi-gas model. The level density 


parameter a is given by 
a(E*, Z, N) = GA)[1 + Ey (Z, N)f(E* - AKE* — A] (9) 


Here, asymptotic Fermi-gas value of the level density param- 
eter 6(4) = aA + B ATP... E (Z, N) is the shell correction 
energy, and f(E”) is the shell damping factor. The parame- 
ters a, f and b, are taken to be 0.114, 0.098 and 1.0, respec- 


tively. The inverse cross section cinv is given by the follow- 


sh 


ing formula : 
aR(1-%), gy, 
ey (€) = e (10) 
0, E<V, 
where R, can be expressed as: 
R, = 1.16[(A — A)? + A15]. (11) 


Here, A, is the mass number of the evaporated particle. For 
proton evaporation, the Coulomb barrier V, is parameterized 
by the following formula : 


V, = [115Z(Z — ZW, + 1.6). (12) 


The survival probability of a superheavy nucleus can be 
expressed as ?l: 


Wr(Eon, x, y, J) m POE qe y, J) 
Z DEJ P DED 
1 Tot Bi, J) IG, J)’ 


(13) 
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where Tiot = Ta + I, + Tp. The fission width I; is given by 
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Fig.2 (color online) Comparisons of calculated evaporation 
residual cross sections with the experimental data in 
40 Ar-?Hf reactions. The solid, dashed, and dash-dotted 
lines represent the calculated evaporation residual cross 
sections in p2n, p3n, and p4n evaporation channels respec- 
tively. The circles, squares, and triangles represent the ex- 
perimental data with error bars in p2n, p3n, and p4n evap- 
oration channels, respectively. The experimental data are 
obtained from Ref. P?! 


the Bohr Wheeler formula, 
1 
T(E*,J) = ———— 
cee 2np(E*, J) 
E'-b-ó- 
| pr(E* — b; — £, J)de 
1 + exp[—22(E* — b; — eK ho] 


(14) 


Here, PUE GO y,J) is the realization probability, 
which is given by the Jackson formulaP""l, E? is the ex- 
citation energy before evaporating the i-th particle. B; is 
the separation energy of the i-th particle. T; is nuclear tem- 
perature before the i-th particle evaporated, obtained from 


T; =  E;la. 


2 Results and discussion 


Firstly, in order to find the compound nuclei with the 
largest cross section in the proton evaporation channel, we 
defined a ratio By B, to evaluate the strength of proton evap- 
oration of the compound nuclei. Fig. 1 shows the calcu- 
lated ratios By B, of nuclei 96-29 Nh, OOO Fl, 283-2968 4. 
and ?96797[ y which denoted by the circles, squares, trian- 
gles, and pentagrams, respectively. It can be seen that the 


ratios B/B, of 285796 Mic are the highest and the fission bar- 
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rier is approximately 2-4 times the proton separation energy, 
which means that proton evaporation is easier for 29-296 Mc 
during the de-excitation process, resulting in a larger proton 
evaporation probability. The ratios B/B, of 286-297F] are the 
smallest, which means that 286-297Fl is hard to evaporate pro- 
tons due to the proton closed-shell at Z = 114. From that, one 
can first synthesize compound nuclei with Z=115, and then 
obtain the Fl isotopes by evaporating one proton and n neu- 


trons. 


To verify the reliability of the DNS model, we calcu- 
lated the ER cross sections in ® Ar-- ? Hf reactions as shown 
in Fig. 2. The p2n, p3n, and p4n evaporation channels 
are denoted by solid, dashed, and dash-dotted lines, respec- 
tively. The experimental data in p2n, p3n, and p4n evapo- 
ration channels are represented by circles, squares, and tri- 
angles, respectively). It can be seen that the calculated 
results are in good agreements with the experimental data, 
but at low incident energy the calculated results are signif- 
icantly lower than the experimental data. This is because 
in the de-excitation process, the evaporation of charged par- 
ticles is based on classical considerations and the quantum 
tunneling effect is not consider. 


In order to produce superheavy nuclei in proton evap- 
oration channel, firstly it is necessary to consider the num- 
ber and order of the evaporated protons. Therefore, in Fig. 
3, we calculated the evaporation width of 480429? Am reac- 
tion in which the Mc compound nuclei can be produced. It 
can be seen that the proton evaporation width is the smallest. 
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Fig.3 (color online) The calculated evaporation width of com- 
pound nucleus ?! Mc produced by *Ca42? Am reaction. 
The solid, dashed, and dash-dotted lines represent the 
widths of fission, neutron, and proton, respectively. 


Compared to neutron evaporation width, the survival proba- 
bility decreases by 1-2 orders of magnitude after each proton 
evaporated. Therefore, in this work, we only considered the 
case of evaporating one proton, and the survival probability 
of multiple protons evaporated is too small, which can be ig- 
nored. Secondly, as the excitation energy decreases, the pro- 
ton evaporation width decreases very quickly, and the ratio 
T/T, decreases. The survival probability will be reduced by 
1-2 orders of magnitude after each proton evaporated. In the 
end, we only consider the case where one proton is evapo- 
rated firstly and then n neutrons are evaporated, and the sur- 
vival probability in other cases is ignored. 


The ER cross sections in *8S+*>Es, °°Ti+?Np, 
34Cr42pa, and ?S4?* Es reactions for producing Fl iso- 
topes are shown in Fig. 4. Solid, dashed, dash-dotted, and 
dotted lines represent pln, p2n, p3n, and p4n evaporation 
channels, respectively. For these reactions, the ER cross sec- 
tions in p2n and p3n channels are relatively larger. One can 
find that the ER cross sections in the reaction SS+ Es are 
the highest among the previous three reactions correspond- 
ing to the same compound nucleus. That is because the mass 
asymmetry for the *°S+*>°Es reaction is the largest, leading 
to the highest fusion probability. The maximum ER cross 
sections of unknown isotopes 7”””'F] are 1.1 pb and 15.1 pb 
at EG, = 38.2 MeV and 39.4 MeV, respectively. 385 beam 
can be produced by the fragmentation of 40Ar at a projec- 
tile fragmentation facility ?9l, and the fusion reactions 38S + 
181 Ta and ?*S + ?95pp with 95S beam were conducted success- 
fully. 


In order to produce more unknown Fl isotopes, the 
ER cross sections in the radioactive beam induced reac- 
tion, ?S4? Es, are also predicted. Four unknown isotopes 
,219^F]. are synthesized with cross sections 3.2, 6.0, 4.0, 
and 0.1 pb at EZ, = 42.2, 28.6, 20.2, and 15.0 MeV, respec- 
tively. One can notice that the ER cross sections in ? S+ 4Es 
reaction in p4n channel are comparable to those in p2n and 
p3n channels unlike the other three reactions. That is be- 
cause the compound nucleus in the reaction 42S+25Es has 
more neutrons than the other three reactions, leading to the 
highest survival probability in p4n channel. If the radioac- 
tive beam facilities are upgraded and the intensity of the ^S 
beam is increased to a high quantity in the future, the reac- 
tion ?S-J* Es can be a promising candidate to approach the 
island of stability. 


Considering the sensitivity of the DNS model to certain 
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Fig.4 (color online) The calculated evaporation residual cross sections in 33S+25Es, 9 Ti429Np, * Cr+ Pa, and * S42 Es reactions. 


Solid, dashed, dash-dotted, and dotted lines represent pln, p2n, p3n, and p4n evaporation channels, respectively. Shadow area 


represents the calculated error. 


parameters, it is necessary to evaluate the inaccuracy of the- 
oretical calculations in order to find out the uncertainty of 
predictions. Among the parameters in the model, the shell 
damping factor Eg plays a vital role in the ER cross sec- 


tions? 


. The shell damping factor cannot be precisely de- 
termined both experimentally and theoretically, and its inac- 
curacy is considered based on the experimental data in the 
Ref. 5l. The uncertainty of the calculation results is shown 
in the shadow in Fig. 4. It can be seen that the uncertain- 
ties of the ER cross sections of Fl isotopes are no better than 
one order of magnitude, and the uncertainties of the optimal 


excitation energy are only in the order of 0.1 MeV. 


3 Summary 


In summary, the production of Fl isotopes in proton 
evaporation channels is investigated within the DNS model. 
Among the isotopes around Z = 114, proton evaporation is 
the easiest for Mc isotopes but the most difficult for Fl iso- 
topes due to the proton closed-shell at Z — 114. Hence, 


one can synthesize the Fl isotopes in the proton evapora- 
tion channel. The calculated ER cross sections by the DNS 
model reproduce the experimental data very well. The syn- 
thesis of unknown isotopes ????!F] is studied via 95S Es, 
SOTi--^5Np, and ?Cr-2?Pa reactions corresponding to the 
same compound nucleus. The ER cross sections of 2?0Fl and 
??!F] in the reaction 38S+255Es are the largest, which are 1.1 
pb and 15.1 pb at EZ, = 38.2 MeV and 39.4 MeV, respec- 
tively. In the radioactive beam induced reaction ?S42*^Es, 
four new isotopes ??!?^F] are synthesized with ER cross 
sections of 3.2, 6.0, 4.0, and 0.1 pb at EZ, = 42.2, 28.6, 
20.2, and 15.0 MeV, respectively. The uncertainties of the 
predictions from the inaccuracy of the shell damping factor 
are calculated, and the uncertainties of the ER cross sections 


of Fl isotopes are no better than one order of magnitude. 
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摘要 : 在 双核 系统 (DNS) 模型 下 系统 地 研究 了 质子 蒸发 道中 产生 超重 核 的 可 行 性 。 由 于 Z=114 aS RENE 
在 , 在 质子 蒸发 道中 合成 Fl 同位 素 是 合适 的 。 在 本 工作 中 我 们 只 考虑 了 先 蒸发 一 个 质子 , 然后 蒸发 nm 个 中 子 的 
情况 , 其 他 情况 由 于 截面 太 小 而 被 忽略 。 与 50Ti+283Np fi Crt Pa 反应 相 比 ,未 知 799291 F] 同位 素 在 9884755 Es 
反应 中 的 产生 截面 最 高 ， 最 大 截面 分 别 为 1.1 和 15.1 pb。 未 来 随 着 放射 性 束 设施 的 升级 ，42S+254Fs 是 接近 稳 
定 岛 的 一 个 有 希望 的 弹 靶 体系 ， 该 反应 中 OLAR 的 产生 截面 估计 分 别 为 3.2、6.0、4.0 和 0.1 pb. 

关键 词 : 双核 系统 模型 ; 超重 核 素 ; 熔 合 反应 ; 质子 蒸发 道 ; 产生 截面 
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